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This talk reviews the recent progress in the extraction of bound-state characteristics from the 
operator-product expansion (OPE) for field-theory correlators, which constitutes the basis of the 
method of QCD sum rules. This progress is mainly related to a deeper understanding of one 
of the key ingredients of the method of sum rules - the effective continuum threshold. The 
latter determines to a great extent the numerical values of the bound-state parameters obtained 
by sum rules. The understanding of properties of the effective continuum threshold allows one 
to formulate a new algorithm for fixng this quantity for various correlators and in this way gain 
control over the systematic uncertainties of the bound-state parameters extracted from sum rules. 
We start with examples from quantum mechanics, where bound-state properties may be calculated 
independently in two ways: exactly, by solving the Schrodinger equation, and approximately, by 
the method of dispersive sum rules. Knowing the exact solution is crucial as it allows us to control 
each step of the sum-rule extraction procedure. On the basis of this analysis, we formulate several 
improvements of the standard procedures adopted in the method of sum rules in QCD. We then 
show the impact of these modifications on the extraction of the decay constants of heavy mesons. 
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1. Introduction 

The method of dispersive sum rules [|I|] is one of the widely used methods for obtaining prop- 
erties of ground-state hadrons in QCD. The method involves two steps: 

(i) One calculates the relevant correlator in QCD at relatively small values of the Eucledian 
time; 

(ii) One applies numerical procedures suggested by quark-hadron duality in order to isolate the 
ground-state contribution from this correlator. These numerical procedures cannot determine 
a single value of the ground-state parameter but should provide the band of values containing 
the true hadron parameter. This band is a systematic, or intrinsic, uncertainty of the method 
of sum rules. 

An unbiased judgement of the reliability of the extraction procedures adopted in the method of sum 
rules may be acquired by applying these procedures to problems where the ground-state parameters 
may be found independently and exactly as soon as the parameters of theory are fixed. Presently, 
only quantum-mechanical potential models provide such a possibility: (i) the bound-state param- 
eters (masses, wave functions, form factors) are known precisely from the Schrodinger equation; 
(ii) direct analogues of the QCD correlators may be calculated exactly. 

Making use of these models, we studied the extraction of ground-state parameters from dif- 
ferent types of correlators: namely, the ground-state decay constant from two-point vacuum-to- 
vacuum correlator [^, the form factor from three-point vacuum-to- vacuum correlator [^], and the 
form factor from vacuum-to-hadron correlator [Q]. We have demonstrated that the standard proce- 
dures adopted in the method of sum rules not always work properly: the true value of the bound- 
state parameter was shown to lie outside the band obtained according to the standard criteria. These 
results gave us a solid ground to claim that also in QCD the actual accuracy of the method may be 
worse than the accuracy expected on the basis of applying the standard criteria. 

We realized that the main origin of these problems of the method originate from an over- 
simplified model for hadron continuum which is described as a perturbative contribution above a 
constant Borel-parameter independent effective continuum threshold. We introduced the notion of 
the exact effective continuum threshold, which corresponds to the true bound-state parameters: in 
potential models the true hadron parameters are known and the exact effective continuum thresh- 
olds for different correlators may be calculated. We have demonstrated that 

(i) The effective continuum threshold is not a universal quantity; it depends on the correlator 
considered (i.e., it is different for two-point and three-point vacuum-to-vacuum correlators); 

(i) The effective continuum threshold depends on the Borel parameter and, for the form-factor 
case, also on the momentum transfer [^, 0, ^, 0]. In QCD, the effective threshold depends 
in addition on the renormalization scale. 

In the recent publications ^ we proposed a new algorithm for extracting the parameters of the 
ground state. The idea formulated in these papers is to relax the standard assumption of a Borel- 
parameter independent effective continuum threshold and to allow for a Borel-parameter dependent 
quantity. This talk explains the details of this procedure and its application to decay constants of 
heavy mesons. 
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2. OPE and sum rule in quantum-mechanical potential model 

Within the standard calculations based on perturbation theory, one start with a free Lagrangian 
and includes interactions order by order. Respectively, the basic object for perturbative calculations 
is the Feynman propagator. However, in a confining theory like QCD, the full propagator in the 
"soft" region of small momenta differs strongly from the Feynman propagator. The Feynman 
propagator of a nonrelativistic particle with the mass m has the form 

Dp(£,P) = - i (2.1) 

- ImE - iO 

Fig. H compares the Feynman and the exact propagators at £ = in a quantum-mechanical potential 
model for the case of the harmonic-oscillator potential V = map-p- jl. This simple model allows 




Figure 1: Feynman vs exact propagator of a confined nonrelativistic particle. 



one to calculate the exact propagator of a confined particle explicitly and illutrates the typical situ- 
ation for a confining potential of a general form. Obviouly, the Feynman and the exact propagators 
are close to each other for large momentum transfers, but differ sizeably in the "soft" region. As 
soon as the "soft" momenta are essential in Feynman diagrams of the perturbation theory, nonper- 
turbative effects in the propagators become essential. At large values of k?-, one obtains 

£>exact(^^) = £>f(^^) + C2|^ + C4|^ + . . . (2.2) 

Nonperturbative effects in the propagator may be described in terms of an expansion in powers of 
OJp' /!?. These nonperturbative effects appear then as power corrections in the correlation functions. 

Let us consider a realistic quantum-mechanical potential model containing both the Coulomb 
and the confining interactions. The corresponding Hamiltonian has the form 

H = Ho + V{r); Ho = ^, V{r) =V,onf{r) - (2.3) 

2m r 

Polarization operator n(£') and its Borel transform Yl{T) [£— )>r, 1/{H — E) ^ exp{—HT), T the 
Borel parameter] are defined via the full Green function G{E) = \/{H — E) [^: 

n(£) = (r/ = 0|^^|r,- = 0), U{T) = {ff = Q\txp{-HT)\fi = Q) . (2.4) 
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The expansion of G{E) and n(£') in powers of the interaction is obtained with the help of the 
Lippmann-Schwinger equation 

G{E) = Go{E) - Go{E)VGo{E) + Gq{E)V Go{E)V Go{E) + . . . , (2.5) 

where Gq{E) = \/{Hq — E). Since the interaction contains now two parts, Vconf{r) and the 
expansion is a double expansion in powers of Vconf and a. E.g., for the case Vconf ('') = 
one easily obtains the corresponding double expansion in powers of a and (oT (see Fig. ^: 



nopE(r) = npert(r)+npower(r) 



npert(r) 



/ m \3/2 
\2%f) 



n 



power 



■T) = i^)"" 



1 + V2KmTa + -mn^Ta^ 

-l«^r^(i + liV2^a)+^co^r^ 



(2.6) 



One can see here a "perturbative contribution" (i.e. the one not containing the confining potential), 
and "power corrections" given in terms of the confining potential (including also mixed terms 
containing contributions from both Coulomb and confining potentials). A perturbative contribution 
may be written in the form of spectral representation [ p^ ] yielding 



noPE(r) = y^/ze-^^Ppert(z)+npower(r), 



m \3/2 



Ppert(^) = (^) 



2J — + V2Kma H —= — 



The "physical" representation for H{T) — in the basis of hadron eigenstates — reads: 




+ 




+ 




+ 



n; 



a- 




(2.7) 



nphys(r) = (r/ = 0|exp(-//r)|r,- = 0) = £/?„exp(-£„r), R„ = \^n{r = 0)\\ (2.8) 

n=0 



Figure 2: Expansion of the polarization functions in powers of the interaction. 
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Sum rule is just the statement that the correlator may be calculated in two ways — using the basis 
of quark states (OPE) or confined bound states — leading to the same result: 

nopE(r) = nphys(r). (2.9) 

In order to isolate the ground-state contribution one needs the information about the excited states. 
A standard Ansatz for the hadron spectral density has the form [|l|] 

Pphys(z) =/?^>5(z-£'g) + 0(z-Zeff)ppert(z). (2.10) 

It assumes that the contribution of the excited states may be described by contributions of diagrams 
of perturbation theory above some effective continuum threshold Zeff- This effective continuum 
threhsold (different from the physical continuum threhsold which is determined by hadron masses) 
is an additional parameter of the method of sum rules. Using Eq. ( |2.10| ) yields 

Zeff 

RgC-'^^'^ = j dze-'^ Pp,ri{?) + npower(r) = Hdual (r,Zeff) • (2.11) 


As soon as one knows Zeff, one immediately obtains estimates for Rg and Eg, Rd\iai{T,Zen) and 
Eduai(T,Zes) = — '^rlogn(r,Zeff)- Thcsc however depend on unphysical parameters T and Zeff- 

3. Exact effective continuum threshold 



Eq. ( 2.10| ) is motivated by quark-hadron duality which claims that far above the threshold the 



hadron spectral density is well described by diagrams of perturbation theory. However, near the 
physical threshold — and this very region turns out to be essentail for the calculation of ground-state 
properties — the duality relation is violated. 



Zexact(T) [GeV] 




— 1 T [GeV"'] 

0.2 0.4 0.6 0.8 1 



Figure 3: Exact effective threshold for the polarization operator in a quantum-mechanical model for the 
potential containing confining harmonic-oscillator and Coulomb interactions (^. 



The advantage of quantum mechanics is that the exact Eg and Rg may be obtained by solving 



Schrodinger equation [|11|]. We can then calculate Zgff from 

Zeff 

Rge-^^^ = J dze-'^PpM+^po^eriT). (3.1) 
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Figures and || provide the exact results for the effective thresholds obtained for the ground 
state of the Hamiltonian (23). We consider a set of confining potentials Vconf('') = (^n {mr)" , n = 
2, 1, 1/2, and adopt parameter values appropriate for hadron physics, i.e., m = 0.175 GeV for the 
reduced constituent light-quark mass and a = 0.3, and adapt the strengths a„ in our confining 
interactions such that the Schrodinger equation yields for each potential the same *F(r = 0) = 0.078 
GeV3/2, ^^^^.^ requires Oi = 0.71 GeV, ai = 0.96 GeV, and ai/2 = 1.4 GeV 0. 

As shown in Fig. |3|, the obtained "exact threshold" Zen{T) for the polarization operator turns 
out to be a slightly rising function of T [^. 



keff(Q) [GeV] 
0.9 



0.8 



Vconf ~ Vconf ~ ^ ^zoni ~ 1"^^^ 




2q x1/3 



kLD=(6;r^Rg) 



2 4 



8 10 



Q [GeV] 



Figure 4: Exact effective threshold ken{T = 0,2) for the elastic form factor for various confining inter- 
actions Vconf('") — <7« (mr)", n — 2,1,1/2 [^] Short-dashed (black) line: harmonic-oscillator confinement, 
n — 2; full (red) line: linear confinement, n = 1; long-dashed (blue) line: square-root confinement, n— 1/2. 
/?g^|*(r = 0)|2. 
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Figure 5: Exact effective threshold feeff (T = 0, Q) for the transition form factor for the case of harmonic- 
oscillator potential [01. Rg = |*(r ^ 0)|2. 



It should be taken into account that in general the effective continuum threhsold depends on the 
type of the correlator considered, on the details of the confining interaction, and for the three-point 
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functions, also on Q^. Figures ^ and || present the results of our calculations of the exact effective 
thresholds for the vacuum three-point correlators at T = 0. Obviously, the exact effective thresholds 
in the region of small momentum transfers Q depend rather sizeably on and on the type of the 
hadron observable under consideration. Interestingly, due to the factorization properties of the 
amplitudes of exclusive processes at large momentum transfers, the exact effective thresholds for 
different correlators approach the same universal constant, related to the value of the bound-state 
wave function at the origin. For further details consult [^. 



4. OPE and sum rule in QCD 

Now, we briefly recall the main ingredients of the QCD sum-rule calculation of the ground- 
state decay constant. The basic object is the two-point function in QCD 

n(/) = ij dxe'P'^ialT (75 W75^(0)) \D.), 75 W = {mQ + m)qiY5Q{x) (4.1) 

The physical QCD vacuum \Q.) is expected to have a complicated nontrivial structure and to differ 



from perturbative QCD vacuum |0). One calculates ( [4.1| ) by constructing the Wilsonian operator 



product expansion (OPE) for T-product of two current operators: 

T(j5{x)jl{0)) =Coix\n)l+Y,Cn{x\n) : d{0,n) : (4.2) 



One now needs to describe properties of the physical QCD vacuum. This is done by introducing 
condensates - nonzero expectation values of gauge-invariant operators over physical vacuum: 

{a\:d{0,n):\Q.)^0. (4.3) 

As the next step, one obtains the dispersion representation for the purely perturbative contribution 
to n(p^) and performs the Borel transform (p^ — )• t) which corresponds to turning from Green 
functions in Minkowski space to time-evolution operator in Euclidean space. The Borel transform 
kills the possible subtraction terms in the dispersion representation for Tl{p^) and improves the 
convergence of the perturbative expansion. Finally, one comes to 

00 

U{t) = J e^'''ppenis,OC,mQ,lJ.)ds + Upoweri'^,mQ,IJ.), 

Ppert(.,M) = p("H^) + ^P^'\s) + f pP)(.) + • • • (4.4) 

n \ n J 

where npower(f)M) is the nonperturbative contribution to n(T). npower(T>M) is given by a power 
expansion in T in terms of the condensates and rad corrections to them. The rum rule reads 

noPE(T) =nhadron(T) (4.5) 

The hadronic part of the sum rule contains the contributions of the ground state and the hadron- 
continuum states. At this point one invokes quark-hadron duality approximation: 

/lOO poo 

/ <isexp(-5T)Ppert(^) = / (i5exp(-ST)Phadr(^)- (4.6) 

■^'^eff '^'^'phys.cont. 
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With the help of this duaUty approximation, the contribution of the excited states cancels against 
the high-energy region of the perturbative contribution, and relates — quite similar to the case of 
quantum mechanics — the ground-state contribution to the low-energy part of Feynman diagrams 



of perturbative QCD and power corrections [ ]12[ ] : 

f^M^e^"Q^= J f'"-"Ppert(^,a,me,;U)t/5 + npower(T,mg,;U) = ndual(T,At,^eff) (4.7) 

{iriQ+m,,)^ 

In order to extract the decay constant one should fix the effective continuum threshold s^n which, 
as illustrated in quantum mechanics, should be a function of t; otherwise the T-dependences of 
the l.h.s. and the r.h.s. of ( p77| ) do not match each other. The exact s^n corresponding to the exact 
hadron decay constant and mass on the l.h.s. is of course not known. The extraction of hadron 
parameters from the sum rule consists therefore in attempting (i) to find a good approximation to 
the exact threshold and (ii) to control the accuracy of this approximation. For further use, we define 
the dual decay constant /^uai and the dual invariant mass M^uai by relations 

/dual(T) =Mg^^''^e^ndual(T,^eff(T)), (4.8) 
MLM = -^logndual(T,Seff(T)). (4.9) 

5. A new algorithm for fixing the effective continuum threshold 

If the mass of the ground state is known, the deviation of the dual mass from the actual mass 
of the ground state gives an indication of the contamination of the dual correlator by excited states. 
Assuming a specific functional form of the effective threshold and requiring the least deviation 



of the dual mass ( |4.8| ) from the known ground-state mass in the T-window leads to a variational 



solution for the effective threshold. As soon as the latter has been fixed, one calculates the decay 



constant from (4.8) 



Our algorithm for extracting ground-state parameters reads: 

(i) Consider a set of Polynomial T-dependent Ansaetze for Seff: 



(ii) Calculate the dual mass for these T-dependent thresholds and minimize the squared differ- 
ence between the "dual" mass Mj^^^ and the known value M| in the T-window 



x'^^X[MiM-Mlf. (5.2) 



!=1 



(n) 

This gives the parameters of the effective continuum thresholds s- . 
(iii) Making use of the obtained thresholds, calculate the decay constant. 
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X n=l 
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Figure 6: The dual mass (left) and the dual decay constant (right) for different Ansatze for the effective 
continuum threshold. 



Figure ^ illustrates the application of our procedure of fixing the effective continuum threshold 
for the realistic case of the D-meson decay constant. Clearly, the T-dependent effective continuum 
threshold leads to a much better stability of the dual mass and allows one to reproduce the knows 
mass of the ground state in a rather broad window of the Borel parameter. 

The crucial question which emerges at this point is how to interprete the results obtained with 
different effective thresholds. To answer this question we compare the extracion of the decay 
cosntant in QCD (where the exact result is not known) with the extraction of the decay constant in 
potential model (where the exact result is known). 

Figure compares the application of our algorithm in quantum mechanics and in QCD [|l3|]. 
In both cases a very similar pattern emerges. In potential model, the exact value of the decay 
constant lies in the band of values provided by the results corresponding to linear, quadratic, and 
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Figure 7: The outcome of our algorithm in potential model (left) and in QCD (right). 
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cubic effective thresholds. We have checked that this result holds independently of the specific 
shape of the confining potential and for a broad range of the parameters of the potential model. 
Since the extraction procedures in potential model and in QCD exhibit precisely the same pattern, 
also in QCD we take this band of values as the estimate of the intrinsic uncertainty of the sum-rule 
estimate for the decay cosntant. So we come to the final point of our algorithm: 

(iv) Take the band of values provided by the results corresponding to linear, quadratic, and 
cubic effective thresholds as the characteristic of the intrinsic uncertainty of the extraction 
procedure. 

Obviously, the T-dependent polynomial effective continuum threshold allows one to reproduce 
the known mass of the ground state much better; this means that the contrinution of the ground 
state from the dual correlator is isolated with a better accuracy. Consequently, the accuracy of the 
extraction of the ground-state decay constant imroves visibly. 

6. Decay constants of charmed mesons 

We now discuss the application of our algorithm to the extraction of the decay cosntants of D 



and Ds mesons [14]. 



6.1 The Borel window 

First, we must fix our working T-window where, on the one hand, the OPE gives a sufficiently 
accurate description of the exact correlator (i.e., higher-order radiative and power corrections are 
small) and, on the other hand, the ground state gives a "sizable" contribution to the correlator. Since 
the radiative corrections to the condensates increase rather fast with T, it is preferable to stay at the 
lowest possible values of t. We shall therefore fix the window by the following criteria Jl^ : (a) In 
the window, power corrections HpowerC'^) should not exceed 30% of the dual correlator HduaK^j'^'o). 
This restricts the upper boundary of the T-window. The ground-state contribution to the correlator 
at this value of T comprises about 50% of the correlator, (b) The lower boundary of the T-window 
is fixed by the requirement that the ground-state contribution does not fall below 10%. 

6.2 Uncertainties in the extracted decay constant 

Clearly, the extracted value of the decay constant is sensitive (i) to the precise values of the 
OPE parameters and (ii) to the prescription for fixing the effective continuum threshold: trying 
different Ansatze for the effective continuum threshold, one obtains different estimates for the 
decay constant. The corresponding errors in the resulting decay constants are called the OPE- 
related error and the systematic error, respectively. Let us discuss these in more detail. 

6.2.1 OPE-related error 

We perform the analysis making use of the MS scheme in which case the OPE exhibits a 
reasonable convergence. The corresponding OPE parameters are summarized here: 

mc{mc) = (1.279±0.013)GeV,m(2GeV) = (3.5 ±0.5)MeV,m,(2 GeV) = (100±10)MeV, 

a5(Mz) = 0.1 176 ±0.0020, <|^GG^ = (0.024 ±0.012) GeV^ 

(^^)(2GeV) = -((267±17)MeV)^(ss)(2GeV)/(^^)(2GeV) =0.8±0.3. (6.1) 
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The value of the OPE-related error is obtained as follows: We perform a bootstrap analysis by 
allowing the OPE parameters to vary over their ranges using 1000 bootstrap events. Gaussian 
distributions for all OPE parameters but pL are employed. For pL we assume a uniform distribution 
in the corresponding range, which we choose to be 1 GeV < /X < 3 GeV for charmed mesons. The 
resulting distribution of the decay constant turns out to be close to Gaussian shape. Therefore, the 
quoted OPE-related error is a Gaussian error. 

6.2.2 Systematic error 

The systematic error of a hadron parameter determined by the method of sum rules (i.e., the 
error related to the limited accuracy of this method) represents perhaps the most subtle point in the 
applications of this method. So far no way to arrive at a rigorous — in the mathematical sense — 
systematic error has been proposed. Therefore, we have to rely on our experience obtained from 
the examples where the exact hadron parameters may be calculated independently of the method 
of dispersive sum rules and then compared with the results of the sum-rule approach. As we have 
seen above, the band of values obtained from linear, quadratic, and cubic Ansatze for the effective 
threshold encompasses the true value of the decay constant; moreover, the extraction procedures in 
quantum mechanics and in QCD are even quantitatively rather similar. Therefore, we believe that 
the half-width of this band may be regarded as realistic estimate for the systematic uncertainty of 
the decay constant. Presently, we do not see other possibilities to obtain a more reliable estimate 
for the systematic error. 



6.3 Decay constant of the D meson 

The T-window for the charmed mesons, T = (0. 1 —0.5) GeV^^, is chosen according to the 
criteria formulated above. Figure |6| shows the application of our procedure of fixing the effective 
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Figure 8: (a) Distribution of fo obtained by the bootstrap analysis of the OPE uncertainties, (b) Summary 
of findings for fo- The estimate obtained with the constant threshold includes the OPE uncertainty only; for 



the T-dependent QCD-SR result the error shown is the sum of the OPE and systematic uncertainties in ( p. 2D 
added in quadrature. 
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continuum threshold and extracting the resulting fj). Figure ^a depicts the result of the bootstrap 
analysis of the OPE uncertainties. The distribution has a Gaussian shape, and therefore the cor- 
responding OPE uncertainty is the Gaussian error. Adding the half-width of the band deduced 
from our T-dependent Ansatze for the effective continuum threshold of degree « = 1,2,3 as the 
(intrinsic) systematic error, we obtain the following result: 



fD = (206.2±7.3(OPE) ±5.1(,y,t))MeV. 



(6.2) 



The main sources of the OPE uncertainty in the extracted fu are its renormalization-scale depen- 
dence and the error of the quark condensate. 

For a T-independent Ansatz for the effective continuum threshold a bootstrap analysis entails 
the substantially lower range /^"~°'' = (l81 .3 ± 7.4(ope)) MeV, which differs from our T-dependent 



result ( |6.2| ) by ~ 10%, i.e., by almost three times the OPE uncertainty. Moreover, as we have shown, 
making use of merely the constant Ansatz for the effective continuum threshold does not allow one 
to probe at all the intrinsic systematic error of the QCD sum rule. From our result ( |6!2[ ) the latter 
turns out to be of the same order as the OPE uncertainty. 

Allowing the threshold to depend on T leads to a clearly visible effect and brings the results 
from QCD sum rules into perfect agreement with recent lattice results and the experimental data 
(Fig. |8|b). This perfect agreement of our result with both experimental data and lattice results 
provides a strong argument in favour of the reliability of our procedure. 



6.4 Decay constant of Ds meson 

We now apply our algorithm to the decay cosntant of D., mesons. Performing the bootstrap 
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Figure 9: (a) Distribution of fo, obtained by the bootstrap analysis of the OPE uncertainties. Gaussian 
distributions for all OPE parameters but /x. For pL we assume a uniform distribution in the range 1 GeV < 
ji <2> GeV. (b) Summary of findings for /o, . Lattice results are from [|l5[ [l^] for two dynamical light flavors 
{Nf = 2) and from [17[ |l8|] for three dynamical flavors {Nf = 3). The triangle represents the experimental 
value from PDG ||l9|. The estimate obtained with the constant threshold includes the OPE uncertainty only; 
for the T-dependent QCD-SR result the error shown is the sum of the OPE and systematic uncertainties in 
(|5.3|), added in quadrature. 
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analysis of the OPE uncertainties, we obtain the following estimate: 

/z), = (245.3 ± 15.7(oPE) ±4.5(,y,t)) MeV. (6.3) 

As in the case of fj), a constant-threshold Ansatz yields a substantially lower value: /^""^ = 
(218.8 ±16.1(oPE))MeV. 

6.5 The ratio /z5,//z) 

For the ratio of the D and decay constants we report the sum-rule prediction 

/d,//d = 1.193 ±0.025(OPE) ±0.007(,y,t). (6.4) 



This value is to be compared with the PDG average foj fo = 1.25 ±0.06 [19] as well as with the 



recent lattice results foj fn = 1.24±0.03 ^ for Nj = 2 and /d,//d = 1.164±0.011 Q and 



fuj Id = 1.20 ±0.02 [|T8p for Nf = 3. The error in (5.4) arises mainly from the uncertainties in the 
quark condensates {ss) / {qq) = 0.8 ± 0.3. 

6.6 Conclusions on decay constants of charmed mesons 

We presented a detailed analysis of the decay constants of charmed heavy mesons with the 
help of QCD sum rules. Particular emphasis was laid on the study of the uncertainties in the 
extracted values of the decay constants: the OPE uncertainty related to the not precisely known 
QCD parameters and the intrinsic uncertainty of the sum-rule method related to a limited accuracy 
of the extraction procedure. 

Our main findings may be summarized as follows. 

(i) The perturbative expansion of the two-point function in terms of the pole mass of the heavy 
quark exhibits no sign of convergence. However, reorganizing this expansion in terms of the cor- 
responding running mass leads to a clear hierarchy of the perturbative contributions. Interestingly, 
the decay constant extracted from the pole-mass OPE proves to be sizeably smaller than the one 
extracted from the running-mass OPE. In spite of this numerical difference, the decay constants 
extracted from these two correlators exhibit perfect stability in the Borel parameter [|T^. This ex- 
ample clearly demonstrates that stability per se does not guarantee the reliability of the sum-rule 
extraction of any bound-state parameter. 

(ii) We have made use of the Borel-parameter-dependent effective threshold for the extraction of 
the decay constants. The T-dependence of the effective threshold emerges quite naturally when one 
attempts to increase the accuracy of the duality approximation. According to our algorithm, one 
should consider different polynomial Ansatze for the effective threshold and fix the coefficients in 
these Ansatze by minimizing the deviation of the dual mass from the known actual meson mass 
in the window. Then, the band of values corresponding to the linear, quadratic, and cubic Ansatze 
reflects the intrinsic uncertainty of the method of sum rules. The efficiency of this criterion has 
been tested before for several examples of quantum-mechanical models. This strategy has now 
been applied to the decay constants of heavy mesons. 
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(iii) We obtained the following sum-rule estimates for the decay constants of the charmed D and 
Ds mesons: 

fo = (206.2±7.3(oPE)±5.1(,yst))MeV, (6.5) 
fo, = (245.3 ± 15.7(oPE) ±4.5(syst)) MeV. (6.6) 

We point out that we provide both the OPE uncertainties and the intrinsic (systematic) uncer- 
tainty of the method of sum rules related to the limited accuracy of the extraction procedure. 
In the case of /o the latter turns out to be of the same order as the OPE uncertainty. Notewor- 
thy, adopting a T-independent effective threshold leads to a substantially lower range f^~^^ = 
(181.3 ib7.4(opE)) MeV, which differs from our T-dependent result ( |63| ) by almost three times the 
OPE uncertainty. 

(iv) Our study of charmed mesons clearly demonstrates that the use of Borel-parameter-dependent 
thresholds leads to two essential improvements: 

a. The actual accuracy of the decay constants extracted from sum rules improves considerably. 

b. Our algorithm yields realistic (although not entirely rigorous) estimates for the systematic 
errors and allows one to reduce their values to the level of a few percent. Due to the application 
of our prescription, the QCD sum-rule results are brought into perfect agreement both with the 
experimental results and with lattice QCD. 

7. Conclusions 

The effective continuum threshold ^eff is an important ingredient of the method of dispersive 
sum rules which determines to a large extent the numerical values of the extracted hadron parame- 
ter. Finding a criterion for fixing Sgft poses a problem in the method of sum rules. 

• i'eff depends on the external kinematical variables (e.g., momentum transfer in sum rules for 3- 
point correlators and light-cone sum rules) and ' 'unphysical" parameters (renormalization scale /i, 
Borel parameter t). Borel-parameter T-dependence of ^eff emerges naturally when trying to make 
quark-hadron duality more accurate. 

• We proposed a new algorithm for fixing T-dependent s^fi, for those problems where the bound- 
state mass is known. We have tested that our algorithm leads to more accurate values of ground- 
state parameters than the "standard" algorithms used in the context of dispersive sum rules before. 
Moreover, our algorithm allows one to probe "intrinsic" uncertainties related to the limited accu- 
racy of the extraction procedure in the method of QCD sum rules. 
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